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Considering the dynamics of a polymer with finite extensibility placed in a chaotic flow 
with large mean shear, we explain how the statistics of polymer extension changes with 
Weissenberg number, Wi, defined as the product of the polymer relaxation time and the 
Lyapunov exponent of the flow. Four regimes, of the Wi number, are identified. One 
below the coil-stretched transition and three above the coil-stretched transition. Specific 
emphasis is given to explaining these regimes in terms of the polymer dynamics. 



Introduction. Recently a number of experimental observations resolving the dynam- 
ics of individual polymers (DNA molecules) in a permanent shear flow have been reported 
by Smith et al. (1999)] see also |Hur et al. (2001)| These ex perimental results and the 



subsequent theoretical/numerical study of |Hur et al. (2000) have focused on the analy- 



sis of the power spectral density and simultaneous PDF of the polymer extension in the 
permanent shear, with fluctuations driven by thermal noise. In another experimental de- 
velopment by Groisman & Steinberg (2000,2001,2004), a chaotic flow state called by the 
authors "elastic turbulence" was observed in dilute polymer solutions. This flow consists 
of regular (shear-like) and chaotic components, the latter being weaker. Resolving an 
individual polymer in this chaotic steady flow was the next challenging but still accessible 



task achieved by Gerashchenko et al. (2004) The coil-stretch transition, predicted by 
Lumley (1973)| (see also Balkovsky et al. (2000) and Chertkov (2000) I, was observed in 



direct single-polymer measurements by Gerashchenko et al. (2004) 



In this letter we present a theoretical analysis of the polymer extension statistics in a 
chaotic flow with large mean shear, s, e.g. of the type corresponding to the elastic turbu- 
lence setup described by Groisman & Steinberg (2000,2001,2004). It is assumed that the 
flow is statistically steady and thus the polymer attains a statistically steady distribution 
as well. We establish the main features of the extension probability distribution, paying 
special attention to the PDF tails. 

The structure of this letter is as follows. We begin introducing the basic dumb-bell-like 
equation governing dynamics of the polymer end-to-end vector, i?, in a non-homogeneous 
flow. Even though the prime interest of this letter is to describe the statistics of the 
polymer extension, R (the absolute value of R) the dynamics of R is tightly linked 
to the angular dynamics. This angular dynamics and the related statistics were the 
subjects of our recent study, see Chertkov et al. (2004) A brief explanation of these 



recent results, relevant to this letter concludes our introduction. We then focus on 
the main subject - to describe the statistics of the polymer extension. We formulate 
the basic stochastic equation governing the dynamics of the polymer extension. Then 
we analyze the structure of the extension PDF which shows a strong dependence on the 
Weissenberg number, Wi. We consider four cases corresponding to qualitatively different 
PDF behaviors. We explain how the typical extension depends on Wi and examine 
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Figure 1. Scheme explaining polymer orientation geometry. 



the tails of the extension PDF, for R less than and larger than typical values. The 
structure of the tails is complicated, consisting in some cases of a number of asymptotic 
sub-intervals. We explain the dynamical origin of all the sub-intervals. To illustrate our 
generic analytical results, we present in Fig. [^graphs, corresponding to the four different 
regimes, obtained by direct numerical simulation made within a model of short-correlated 
velocity statistics and of the so-called FENE-P modeling of polymer elasticity. 

Model. We consider a single polymer molecule which is advected by a chaotic / turbulent 
flow (i.e. the polymer moves along a Lagrangian trajectory of the flow) and is stretched 
by velocity inhomogeneity. The polymer stretching is characterized by the molecule's 
end-to-end separation vector, i2, satisfying the following dumb- bell- like equation (see 
e.g. 



Hinch (1977) Bird et al. (1987) I 



^tR^ = RJV,v,--^{R)R, + C^■ (0.1) 

Here 7 is the polymer relaxation rate and Qi is the Langevin force. The velocity gradient 
V jVi is taken at the molecule position. The velocity difference between the polymer end 
points is approximated in Eq. H0.1|l by the first term of its Taylor expansion in the end- 
to-end vector. This approximation is justified if the polymer size is less than the velocity 
correlation length. The relaxation rate 7 in Eq. (|0.1() is a function of the extension R 
which varies from zero upto a maximum value iimax corresponding to a fully stretched 
polymer. We assume that the relaxation is Hookean for R <C Rma^, i-e. 7(^) is well 
approximated by a constant 70 there, while it diverges (the polymer becomes stiff) for 

R ^ Rma,K • 

We focus on the situation in which the effect of velocity fluctuations is stronger than 
that of thermal fluctuations, so that the Langevin force C in Eq- l|n-l|l can be neglected. 
We consider the case in which the steady shear flow is accompanied by weaker random 
velocity fluctuations. This is also the setting realized in the elastic turbulence experiments 
by Groisman & Steinberg (2000,2001,2004). We choose the coordinate frame associated 
with the shear flow, as shown in Fig. ^ where the mean shear flow is characterized by the 
velocity (sy, 0, 0) and s is positive. Then the polymer end-to-end vector R is conveniently 
parameterized by the spherical angles and 9: Rx = R cos 9 cos (j), Ry = Rcos9 smcj), 
Rz = Rsm9. In terms of these variables, Eq. 1)0. l|l (with the Langevin term omitted) 
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transforms into the foUowing set of equations: 

dt(t) = -ssm^(j) + ^^, (0.2) 
dtO = — s sin (p cos (f> sin 9 cos 6 + £^0 , (0-3) 
dt\nR = -j{R) + s cos^ 9 cos sin + C|| , (0-4) 

where £s and ^|| are random variables related to the fluctuating component of the 
velocity gradient. Note, that the angular (orientational) dynamics described by Eqs. 
HU.2IU.3|I decouples from the dynamics of the polymer extension, R. Another remark is 
that, at 7 = 0, Eq. (|0.4|) describes the divergence of neighboring Lagrangian trajectories. 
The character of the angular dynamics (closely related to the Lagrangian dynamics 



in the flow) was discussed in detail by Chertkov et al. (2004) Typically, the polymer 



orientation fluctuates near its preferred direction determined by the shear. Characteristic 
values of the fluctuations, both in 9 and 0, are determined by the average value of the 
angle 0, (j)f , while the average value of 9 is zero. (We consider parametrization with 
the angles taken inside the torus —tt/2 < (j),9 < 7r/2.) The value of is small due to 
assumed weakness of the velocity gradient fluctuations ^ in comparison with the shear rate 
s. Note that (j>t is directly related to the value of the Lyapunov exponent of the flow. A: 
A = s(f)t ■ We make the natural assumption that the flow velocity is correlated at the time 
Tt = X~^. The random terms ^0 and £e in Eqs. (|().2I().3|I are relevant (i.e. comparable to 
their deterministic counter-parts) only in the narrow angular region \9\ < (pt. Outside 
this "stochastic domain" the angular dynamics is mainly deterministic, i.e. the effect of 
the stochastic terms, and ^g, is not essential there. The deterministic motion leads 
to polymer flipping (i.e. reversing its orientation). The flips interrupt slower stochastic 
wandering near the shear-preferred direction. The process of the deterministic/stochastic 
regimes alteration (tumbling) is a-periodic. 

The Statistics of Polymer Extension. We consider the case of a statistically steady 
random flow, that leads to stationary statistics as a function of the angles, 9, </>, and of the 
extension R. Our goal here is to describe the statistics of R that emerges as a result of 
the balance between elasticity-driven contraction and extension, caused by fluctuations 
in the flow. 

For the principal part of the dynamics, the basic dynamical equation Eq. 1)0. 4|l can 
be simplified. First, the main contribution to the R dynamics stems from the region 
of small angles where sin0 can be replaced by (j) and cos^ by unity. Second, the term 
^11 is potentially important only in the stochastic region where competes with s^'^. 
However, assuming ~ ^||, we conclude that ^|| is negligible in comparison with scj) 
there. Therefore one arrives at: 

dt\nR^--^{R) + s(l). (0.5) 

Note that Eq. (|0.5|) is inapplicable when R is close to its minimum value during a flip 
(since the angles 0,0 are of order unity there). Another remark is that Eq. (|0.5|) is 
correct for R Rt {Rt is the typical length of the polymer in the absence of the flow) 
where the Langevin force can be neglected. 

The statistics of R is determined by the interplay of the two terms on the right-hand 
side of Eq. H0.5|) . Since the average value of s(j) is equal to the Lyapunov exponent, 
A, the dimensionless parameter characterizing the statistics of the polymer extension 
is the Weissenberg number, Wi = A/70, which grows with the strength of the shear, 
or/and, of the velocity fluctuations. At Wi = 1, when the two terms on the right hand 
side of Eq. 10. 5|) balance each other, the system undergoes the so-called coil-stretched 
transition, sec Lumlcy (1969)1 |Lumley"[T97"3),Balkovsky et al. (2000)| [Chertkov (2000) 
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Figure 2. PDF of the polymer extension, R, measured in the units of maximal extension, for 
different values of the Weissenberg number, Wi, obtained from numerical simulations of the 
stochastic equations explained in the text. 



and Balkovsky et al. (2001) for details. We find, however, that in the specific case of 
strong shear considered in the letter additional qualitative changes in the PDF of R 
occur at > 1 so that the overall picture is richer than the case of isotropic velocity 
statistics. Below we describe characteristic features of the extension PDF as a function 
of Wi. 

To illustrate our generic analytical results we plot in Fig. |21four graphs of the extension 
PDF obtained by numerical simulations based on modification of Eq. H0.5|l . dt\nR = 
—7 + ssin0cos0, which allows the correct reproduction of the flips, and Eq. H0.2I) with 
the stochastic term chosen to be (5-correlated in time. The simulations were done with 
7(i?) = 7o/(l — i?^/i?^a^x)> corresponding to the so-called FENE-P model of the polymer 



elasticity, see e.g. Bird et al. (1987) 

Wi < 1, a > 0. We begin by discussing the case Wi < 1, for which the polymer is only 
weakly stretched. Then, typically, the molecule stays in the "coil" state characterized 
by the thermal size, Rt, which emerges as the result of a balance between the Langevin 
driven extension and contraction (relaxation) related to polymer elasticity. Let us recall 
that due to a large number of monomers in the polymer molecule, the thermal noise 
induced length, Rt, is much smaller than the maximal polymer extension, i?max- 

At the scales larger than Rt the thermal noise is irrelevant and the extension dynamics 
is described by Eq. (|0.5|l . Being interested in the statistics for large, R » Rt, for 
deviations from the typical value, Rt, one comes to a problem that was already analyzed 
in detail by [Balkovsky et al. (2000)1 [Chertkov (2000)| [Balkovsky et al. (2001)| in which 
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it was shown that the extension PDF has an algebraic tail: 

P{R) oc R-^-" , (0.6) 

with a > 0. Eq. (|0.t)|l holds at i?niax ^ i? 3> Rt where 7(i?) weakly deviates from 
7o. The situation is reflected in Fig. where the algebraic tail is clearly seen. The 
positive value of the exponent a in Eq. I|U.6|I guarantees that the normalization integral 
J dRP{R) converges in the region R ^ Rt- Thus, the normalization coefficient in Eq. 
IjU.BI) is ~ i?^". The exponent a decreases as the Weissenberg number, Wi, increases, 
and it crosses zero at the coil-stretch transition, where Wi = 1. 

The algebraic tail (|0.6() is related to a long (compared to the correlation time Tt) 
process of polymer extension from the typical value Rt to the current value of the 
extension, R ^ Rt- Note that this long extension docs not mean that (during the 
process of extension) the right hand side of Eq. (|0.5(l is always positive since (p fluctuates 
and s4> is larger than 7 only on average. Moreover, the process consists of alternating 
stochastic and deterministic portions (polymer flips during the later ones). In spite 
of the fact that R decreases during the first half of the flip, the initial extension is 
restored (returns to its initial value) during the second half of the flip. Overall, the 
flips do not influence the extension process. The probability W of this a- typically long 
extension process depends exponentially on its duration T, since is a product of 
independent probabilities, each representing a sub-processes of duration Tf. This gives 
the following estimate: Iniy ~ —T/ti. On the other hand, in accordance with Eq. (|0.5() . 
ln(i?/i?T) T/rt- Combining the two estimates, we arrive at the algebraic tail (|().6|l for 
the PDF of extension, P = dW/dR. 

Wi > \, —1 < a < {). Above the coil-stretch transition, when A exceeds 70, the poly- 
mers become strongly extended. In this stretched state the typical size of the polymer, 
i?*, is much larger than Rt- Considering the stationary average of Eq. H0.5|) one finds 
that 7(i?,) = A. 

The left tail of the PDF, corresponding to Rt R ^ R*, is governed by the same 
algebraic law 1)0. However, now a < 0, meaning that the normalization integral 
/ dRP{R) gets a major contribution for R ^ R^,. Therefore, restoring normalization, 
one derives P ~ for the i? ^ i?* tail. The transition from positive to negative 

a corresponds to important change in the nature of the dynamical configuration corre- 
sponding to the algebraic tail: extension, as a typical process for a > 0, is replaced by 
contraction for negative a, so that initially typical extension, ^ i?*, contracts through a 
long, T ^ Tt, (multi-tumbling) evolution. (Physical arguments, clarifying the origin of 
the algebraic tail, are identical to the ones presented above for Wi < 1.) 

The right tail, corresponding to extreme extensions, i?max — R ^ R — R^,, can also 
be explained in some general terms. The domain of extreme deviation is characterized 
by extremely fast relaxation, so that the term on the left hand side of Eq. 1)0. 5|l can be 
neglected. As a result, R and (j) are related to each other locally: 7(i?) = scj). Moreover, 
one finds that, because of the fast nature of the polymer relaxation in the extreme case, 
(p simply follows the respective random term in Eq. H0.2|l . i.e. the term on the left 
hand side of Eq. (|0.2|l can also be neglected resulting in, s(j>^ = C^- other words, the 
extreme configuration is produced through fast anti-clockwise revolution of the polymer 
to a large (in comparison with the typical value (pt) positive angle, 1 3> S> 0f, driven 
by anomalous fiuctuations in Recalculating the PDF of into P{R) one arrives at 

P{R) = 2s-'^j'P^{jys), (0.7) 

where P^ is the simultaneous PDF of the velocity gradient term Note that the 
asymptotic expression given by Eq. H().7|l is not restricted to the case considered in this 
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subsection but applies generically to the description of extreme polymer extensions in all 
regimes. 

Wi > 1, a < — 1. Once a crosses —1, i?* becomes maximum of the extension PDF, 
P(i?). This modification in the PDF shape is accompanied by the emergence of a plateau 
on the left from the maximum (see Fig. 13;), associated with an additional contribution 
to the PDF related to the deterministic angular dynamics. 

Let us explain the origin of the plateau. For angles which are smaller than unity but 
larger than (j>t, R and </> satisfy, 9ilni? = s(/) and dt(l> = ~s(f>'^, as follows from Eqs. 
H0.2I0.5|I . Integrating these equations one arrives at, R = A\t — to] where to and A 
are constants, the latter being estimated by A ~ R^/Tf. Taking into account the fact 
that the time Iq is homogeneously distributed (due to the assumed homogeneity of the 
velocity statistics), and recalculating the measure dt^ into the PDF of _R, one arrives 
at P{R) = C /R^, (C is an i?-independent constant of order unity) corresponding to the 
plateau seen in Fig. 

The "deterministic" contribution to P{R), ~ ^/R*, discussed above, does not cancel 
the "stochastic" one ~ R"R~^~°' corresponding to the long contraction which starts at 
i?,. In fact these contributions co-exist. One finds that for —1 < a < 0, the stochastic 
contribution dominates, in full agreement with the above discussion of Fig. 13). When 
a becomes smaller than —1, the situation is reversed and the deterministic contribution 
dominates. 

The plateau extends from i? i?* down to i? R*4't, where i?*0t is the smallest value 
of R one can achieve under the condition that when the flip begins the initial extension 
is i?*. Note however, that if the initial extension is smaller than i?*, a deterministic 
flip could bring it to values that are smaller than i?*(/)t. Therefore, to understand even 
smaller values oi R, R < R^,(j)t, one should consider flips which begin with an anomalously 
small initial value Rq, Rq < i?» (prepared by some preliminary and long stochastic 
processes, of the type discussed above). The probability density of achieving i?o during 
the preparatory stage is estimated according to Eq. (|0.6|l : ~ R"Rq^~°'. On the other 
hand, any Rq that lies between R and R/4>t transforms dynamically as a result of a fast 
flip into the current value of extension R with the same i?-independent probability, which 
now can be estimated as ~ 1/Ro- Therefore, one arrives at the following estimate for 
the PDF of i? at i? < E*0t: 



Eq. (|0.8I) explains the probability decrease at the smallest R seen in Fig. 

Let us observe the bump in Fig. 13; separating the region of the plateau and the 
smallest R region of the probability decay. To explain the bump, one simply needs to 
account for the angular nonlinearity (with respect to (f> and 9) in the estimate of the 
plateau just discussed. 

Wi :S> 1, a -^i —1. The larger Wi is, the closer i?* approaches i?max- Then the condi- 
tion of fast relaxation, Ry{R) ^ 7, which has already allowed us to analyze the extreme 
asymptotics 1)0. 7|l . also applies to the region in the vicinity of i?, . The smallness of the 
ratio 7/(i?7') suggests that the left hand side of Eq. (|().5|l can be replaced by zero. 
One arrives at 7(i?) — scf), which makes it easy to express the PDF of R through the 
simultaneous PDF of 0: 



where it is also assumed that ^(R)/s < (pt- In this special domain of R and 0, variations 
of Pij, are slow, so that the main dependence on R in Eq. H().9|l is due to the factor y{R). 




(0.8) 



P{R)^s~'j'{R)P4j/s), 



(0.9) 
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Eq. (|().9|l applies to the left (for smaller values of R) of i?* provided the parameter 
i?7'(i?)/7 is large. In this domain, the PDF can be estimated as P ~ ^'{R)/^{R^). At 
even smaller values of i?, the PDF has a plateau, P ^ 7o/[-R*7(i?*)] (which generalizes 
our previous result to the Wi 1 case). This explains the complex behavior of the PDF 
of R shown in Fig. 

Conclusions. In this second work devoted to the statistics of polymer molecules in 
a chaotic flow with mean shear, we focused on analyzing the statistics of the polymer 
extension R. The PDF of R demonstrates complex and rich W^i-dependent behavior 
related to a number of distinct processes governing polymer dynamics. We observed that 
the typical value of the extension associated with the stochastic wandering of the polymer 
orientation around the special shear-preferred direction increases with Wi. In the region 
of maximal stretching, near -Rmax, the major contribution to the PDF originates from 
processes characterized by fast adjustment of the polymer extension to the current value 
of the polymer's angular degree of freedom. We also identified special contributions to 
the PDF tails associated with fast (deterministic) flips and long (stochastic) extension 
or contraction processes. Encouraged by qualitative agreement of our results with the 
newest experimental data of Gerashchenko et al. (2004]} we anticipate that the rich zoo 
of theoretical predictions presented in the paper will be helpful for guiding and testing 
future experimental work in this field. 
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16520a (IK,VL and KT), and a grant from the Dynasty Foundation (KT) are acknowl- 
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